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Based on the R-parity violation option of the minimal supersymmetric Standard Model, we 
examine the scenario where the massive gravitino, relic from the hot big-bang, is the lightest 
supersymmetric particle and can decay through one or several of the trilinear R-parity violating 
interactions. We calculate the rates of the gravitino decay via the various three-body decay 
channels with final states involving three quarks and/or leptons. By taking into account 
the present constraints on the trilinear R-parity violating coupling constants and assuming 
the gravitino and scalar superpartner masses do not exceed ~ 80TeV, it turns out that the 
gravitinos could easily have decayed before the present epoch but not earlier than the big-bang 
nucleosynthesis one. Therefore, the considered scenario would upset the standard big-bang 
nucleosynthesis and we conclude that it does not seem to constitute a natural solution for the 
cosmological gravitino problem. 

1 Introduction 

In supergravity theories 0, the gravitino, namely the spin-3/2 supersymmetric partner of the 
graviton, weakly interacts with all the particle species (including itself) due to the small grav- 
itational strength coupling \JGn = 1/Mp, Gn being the gravitational constant and Mp the 
Planck scale. Hence, the gravitino-gravitino annihilation rate is extremely small so that the 
gravitinos should decouple at an early epoch of the universe history, and moreover at an epoch 
characterized typically by a temperature higher than the gravitino mass: kT^ > m 3/ / 2 ik 
being the Boltzmann constant) 0. Therefore, the relic abundance of the gravitino should be 
large, which is often denoted in the literature as the cosmological "gravitino problem" . One of 
the first solutions for the gravitino problem to be envisaged is to compensate the large gravitino 
relic abundance by a gravitino mass sufficiently small, namely 7713/2 < l&eVu, to respect the 

limit on the present universe energy density: Qq ~ 1. For heavier gravitinos, a second type 
of available solution is by shortening their lifetime so that they do not survive out at the late 
epochs 0. This can be realized in two characteristic options: either the gravitino is not the 
Lightest Supersymmetric Particle (LSP) and thus can decay into an odd number of superpart- 
ners through gravitational and gauge interactions (both of these ones couple an even number of 
superpartners), or it is the LSP. A gravitino LSP can be realized within various supersymmetric 
models, including the gauge mediated supersymmetry breaking models a>B, the models of low 
fundamental energy scalei and even the conventional gravity mediated supersymmetry breaking 
models (for some specific set of the supersymmetry breaking parameters). If the gravitino is the 
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LSP, it can decay only into the ordinary particles of the Standard Model. Such a decay channel 
must involve both gravitational and the so-called R-parity symmetry l'i violating interactions 
(the latter ones couple an odd number of superpartners) . The R-parity violating (ftp) interac- 
tions are written in the following superpotential, in terms of the left-handed superfields for the 
leptons (L), quarks (Q) and Higgs of hypercharge 1/2 (H) and the right-handed superfields for 
the charged leptons (E c ), up and down type quarks (U C ,D C ), 

W ft p = ]L \ 7^ijk L i L j E k + Kjk L iQj D k + 2Kjk U iD C j D k + HiHL^j , (1) 

^i5ik 

i,j, k being flavor indices, A^, A- jfc , A" jfc dimensionless coupling constants and \x\ dimension one 
parameters. Note that in the scenario in which the gravitino is not the LSP, the gravitino 
preferentially decays into an ordinary Standard Model particle and its superpartner through 
gravitational interactions, since the ftp coupling constants are severely constrained by the low- 
energy experimental bounds obtained at colliders 00. 

In the scenario of an unstable gravitino heavier than the LSP, the gravitino decay can 
produce an unacceptable amount of LSP which conflicts with the observations of the present 
mass density of the universe tL3. The scenario containing an unstable LSP gravitino, having 
decay channels which involve ftp coupling constants, is based on the violation of the R-parity 
symmetry. Now, neither the grand unified theories, the string theories nor the study of the 
discrete gauge symmetries give a strong theoretical argument in favor of the conservation of 
the R-parity symmetry in the super symmetric extension of the Standard ModelH3. Hence, the 
scenario with an unstable LSP gravitino constitutes an attractive possibility which must be 
considered as an original potential solution with respect to the cosmological gravitino problem. 

Our main purpose in the present work is to determine whether the scenario of an un- 
stable LSP gravitino decaying via ft p interactions constitutes effectively a natural solution to 
the cosmological gravitino problem. We will concentrate on the trilinear ft p interactions, namely 
XijkLiLjEf,, X' i j k LiQjDl and X'^UfD^D^. With the goal of enhancing the gravitino instability, 
we will consider an optimistic type of scenario in which several trilinear ft p coupling constants 
have simultaneously non-vanishing values. We will not consider the bilinear ft p term jjiiHLi 
of Eq.([j]) which can be rotated away by a suitable redefinition of superfields, for some specific 
choices of the ft p soft supersymmetry breaking termsl 10 H 13 l. The bilinear ft p interactions as well 
as the possible alternative of a spontaneous breaking of the R-parity symmetry have been con- 
sidered, within the context of the cosmological gravitino problem, in a recent study which 
examines the twodpody gravitino decay mode into photon and neutrino. We also mention here 
the workstL3 andli^l in which have been considered the ft p decay reactions of single nucleon into 
a gravitino and a strange meson. 

2 Constraints in a scenario containing an unstable gravitino 

The equilibrium condition, 

r > H(T), (2) 

r being the gravitino decay rate and H(T) the expansion rate of the universe, is reached at an 
epoch at which the cosmic energy density is dominated by the gravitino energy density Pz/2°- 
Therefore, in determining the temperature T 3 / 2 at which the equilibrium condition of Eq.(|2|) is 
reached (the gravitino decay temperature) by solving the equation T = H(T 3 / 2 ), on e must take 
the following expression for the expansion rate, 



H(T) PS \J^Y^p 3 /2, 
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Figure 1: Feynman diagrams of the gravitino three-body decay processes involving the XijkLiLjE% interactions. 
These interactions are represented by the circled vertex. The other vertex correspond to gravitational couplings. 
The arrows denote the flow of momentum of the associated particles. Finally, e, e, v, v and G represent respectively 
the charged leptons, the charged sleptons, the neutrinos, the sneutrinos and the gravitino, i,j, k are generation 
indices, the L, R indices stand for Left, Right (chirality) and the c exponent indicates a charge conjugated particle. 
We have not drawn the charge conjugated reactions which can be trivially deduced from these ones. 



where p 3/ / 2 = (3C(3) /n 2 )(g(T) / ' g(T,i))m 3 /2(kT) 3 , Q(x) being the Riemann zeta function (C(3) = 
1.20206...) and g(T) counting the total number of effectively massless degrees of freedom. 
Hence, T 3 / 2 reads as, 

vr qlTA \V 3 /T 2 Ml\V 3 



kT ^<mM}) • (4) 

At the temperature T 3 / 2 , most of the gravitinos decay. After the gravitinos have decayed and 
their decay energy has been thermalized, the temperature rises to the value kT^, 2 given by, 

, T , „ (t»0<(3))'/ 4 ( mwQtfyrf y/* ( 45 \ l " g{T d ) l l* 

In order to avoid large relic abundance of an unstable gravitino, this one must decay before 
the present epoch. This requirement imposes the following bound on the increased gravitino 
decay temperature TL 2 , 

T' zj2 > 2.75K (kT^ /2 > 2.36 W~ 10 MeV). (6) 

Furthermore, the decay of the gravitinos causes an increase of the temperature, which leads 
to an increase of the entropy density s = (2ir 2 /A5)g(T)(kT) 3 and hence to a decrease of the 
baryon-to-entropy ratio B = tib/s, ub being the baryon number density. Therefore, if the 
gravitinos decay after the nucleosynthesis epoch, the baryon-to-entropy ratio during the nucle- 
osynthesis epoch must be much greater than the present one, leading to the production through 
nucleosynthesis of too much helium and too little deuterium (compared with the constraints on 
the primordial abundances derived from observational data) □. In conclusion, if the gravitinos 
decay before the present epoch, the temperature Ty 2 , reached after the gravitino decay energy 
has been thermalized, has to be higher than the nucleosynthesis temperature, namely, 

kT' zj2 > OAMeV. (7) 

3 Gravitino life time 

The gravitino decay processes which involve the trilinear ft p interactions of Eq.(||) and have the 
dominant rates are obviously the decays through a virtual scalar superpartner into three ordinary 
fermions involving one gravitational and one fi p coupling. Those processes are represented in 
Fig.|l| for the R p interactions of type XijkLiLjE'f,. The reactions involving the A^- fc and X'^ k 
couplings have similar structures. 



As a function of the relevant ratio of squared masses of the gravitino and exchanged sfermions, 
z = my 2 /rh 2 , we find that the partial gravitino lifetime, associated to the decay into a given 
flavor configuration of three light fermions final state, decreases monotonically in the interval 
< z < 1 as, 

- * 109 - 1011 (i)(^) 3 - c < « 8) 

where A stands for any one of the trilinear coupling constants. For light gravitinos of mass 
m 3/2 = 0(GeV) and A^ = 0(1), we see from Eq.(|8|) that the gravitino lifetime exceeds by a 
few order of magnitudes the age of the universe to — 3.2 10 17 sec. 

4 Numerical results and discussion 

In Fig]2| (plain line), we show the values of the gravitino and superpartner masses for which 
kTy 2 > OAMeV. These values have been obtained by assuming that the single dominant 
ftp coupling constant is A 213 and is equal to its present limit 0. 

In Fig.[| we have set all the scalar superpartner masses to a common value, denoted by rh. 
Moreover, the combination of statistical factors ^(Ti) 1 ' 4 /^^^) 1 ' 2 , which enters Eq.(|5|), has 
been set to unity \ Finally, the gravitino decay rate has been multiplied by a factor of 2 in 
order to count the charge conjugated gravitino decay process and the masses of the final state 
particles have been taken into account in the computation. 

We see in Fig.|2| (plain line) that as the superpartner mass increases, larger gravitino masses are 
needed to have kT^, 2 > OAMeV. The reason is that T^ 2 increases with the gravitino mass but 
is suppressed if the superpartner is getting heavier (since Ty 2 is proportional to the squared 
root of the gravitino decay rate, as shown in Eq.@). 

Let us consider now a type of scenario in which several fl p coupling constants, having some 
of the weakest present bounds, are simultaneously non-vanishing and are equal to their present 
limit obtained in the single dominant coupling hypothesis. Based on the strongest constraints 
on the products of fi p coupling constants and on the reviewlll of the present limits on the single 
ftp coupling constants, we find that the scenario of this type leading to the highest gravitino 
decay temperature TL 2 (for m 3 / 2 = ITeV and rh = 1.5TeV) corresponds to the case where the 
simultaneously dominant ft p coupling constants are, 

^132) ^211 J ^223 > ^311 ' ^121 an d A233, (9) 

which have the present bounds (obtained in the single dominant coupling hypothesis) A^ 32 < 0.34 
for TTiq = WOGeV (and for instance A' 132 ~ 1.2 for m~ q = ITeV) 00, A' 211 < 0.06(mj /lOOGeV) 
Bill, A 223 < 0.18(m^ fl /100GeV)il, A?m < 0.10(m^ /lOOGeV)0§l, Am < 0.05(m e ~*/100GeV) 
00 and A 233 < O.O6(m ffl /lOOGeV)0Ilf We note that no \'! jk coupling can be added to the 
set ^ of simultaneously dominant fi p couplings, since the experimental constraints on the proton 
decay rate force any product A^A"^/ to be smaller than 10 -9 , in a conservative way and for 
squark masses below ITeV. This result has been obtained in§l by calculating the proton decay 
rate at one loop level. In contrast, we have checked that no strong constraints exist on any 

Let us assume tentatively that the thermalized degrees of freedom at the gravitino decoupling and decay 
temperatures are the same as those for the minimal supersymmetric Standard Model and the present epochs, 
respectively. Using then the values g(T d ) ~ 915/4 ~ 228.75 and g(T 3/2 ) = 43/11 = 3.909, one obtains: T 3/2 oc 
ffC^d) /ffCfs/a) — 1-96. The actual numerical value of this factor might be larger but not by very much. 

c In Fig.g, the ]ft, p coupling constants of set [] have been set to their perturbativity limit, obtained from the 
requirement of perturbativity up to the gauge group unification scale Ej, since, in the whole intervale of rh 
covered by the figure, this latter limit is more severe than the corresponding present bound, which is a low-energy 
experimental constraint. 
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Figure 2: Domains of the m i f 2 {GeV / e 2 )-rh(GeV / c 2 ) plane (gravitino versus superpartner mass) in which the 
gravitino decay temperature is higher than the nucleosynthesis one, namely kT^ 2 > OAMeV . The region situated 
above the plain line corresponds to kTy 2 > OAMeV in case the dominant lfi v coupling constant is A213 = 1.25. 
The domain situated above the dashed line corresponds to kTy 2 > OAMeV in case the dominant ft p coupling 
constants are \' 132 = 1.04, X' 211 = 1, X' 223 = 1.12, A 311 = 1.12, A121 = 1 and A233 = 1. Finally, the colored region 
corresponds to the situation m 3 / 2 > m which must be considered within a scenario where the gravitino is not the 

LSP. 

product taken among the fl p couplings of the set ||. 

In this optimistic situation where several ft p coupling constants are simultaneously present, the 
gravitino masses corresponding to kT^ 2 > OAMeV are smaller than in the considered case of 
a single dominant fl p coupling constant (see Fig.^), since Ty 2 is typically proportional to the 
fi p coupling constants and decreases with the gravitino mass. 

The conclusion about Fig.^ is that the gravitino mass, and thus the scalar superpartner 
masses, must exceed values of O(80TeV), even within the most optimistic scenarios, so that the 
condition of Eq.(|7]) can be fulfilled. Nevertheless, this requirement conflicts with the constraint, 
fh ~ 0(TeV), coming from the "hierarchy problem", namely the problem of natural coexistence 
of the electroweak symmetry breaking scale and the scale of new physics underlying the Standard 
Model (grand unification scale, string scale . . . ). 

5 Conclusion 

Along with the existence of the cosmic microwave background, big-bang nucleosynthesis is one 
of the most important predictions of the big-bang cosmology. Furthermore, if one assumes that 
the light nuclei (atomic number less than 7) have effectively been produced through the big- 
bang nucleosynthesis, one finds that the theoretical predictions on the abundances of these light 
nuclei are in good agreements with the observational data El Now, as we have seen above, 
if one believes that the light elements have been synthesized through the standard big-bang 
nucleosynthesis, the temperature Ty 2 > reached after the thermalization of the gravitino decay 
energy, must be higher than the nucleosynthesis temperature. Therefore, since we have found 
that this cannot happen in the scenario characterized by an unstable LSP gravitino having a 
decay channel which involves trilinear fi p interactions, this scenario does not seem to provide a 
realistic solution to the large relic abundance of the gravitino. 
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